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An  algorithm  is  derived  Car  determining  for  the  location  of  a  stationary  emitter  based  upon  a 
number  of  angular  measurements  from  a  moving  platform.  It  is  referred  to  as  ihe  asymptotic  algonthar 
and  is  particularly  useful  when  there  is  appreciable  error  in  the  angle  measurements.  It  l»  restricted 
to  platform  motion  which  generates  a  base  tine  small  in  comparison  to  the  range  from  the  platform 
to  the  emitter.  Under  there  conditions  better  retails  arc  obtained  than  from  the  slaadajd  least 
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squares  distance  algorithm.  This  is  accomplished  with  virtually  the  same  amount  of  memory  and 
computer  logic. 

An  algorithm  that  obtains  emitter  velocity  as  well  as  position  is  also  given  in  the  report.  This 
algorithm,  which  requires  a  nonlinear  maneuver  by  the  platform,  gives  useful  results  when  the 
system  has  good  direction-finding  accuracy  and  a  capability  for  recording  the  time  of  a  measurement. 
Typical  results  for  both  algorithms  are  presented  in  tables  from  which  it  is  possible  to  judge  the  merits 
of  the  method  for  a  particular  application. 
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EMITTER  LOCATION  USING  BEARING 
MEASUREMENTS  FROM  A  MOVING  PLATFORM 

INTRODUCTION 

This  report  presents  a  set  of  related  mathematical  algorithms  for  determining  the  location  of  an 
emitter  based  upon  a  number  of  angular  measurements  from  a  moving  platform.  For  all  except  the  last 
algorithm  given,  the  emitter  is  stationary  with  respect  to  the  bearing  measurement  platform.  These 
algorithms  are  particularly  useful  when  there  is  appreciable  error  in  the  angle  measurements,  but  they 
are  restricted  to  platform  motion  which  generates  a  base  line  small  in  comparison  to  the  range  from  the 
platform  to  the  emitter.  An  alternative  approach  to  this  problem  is  given  in  a  paper  by  Mahapatra  (1). 

Before  presenting  the  new  algorithms,  another  will  be  outlined  for  purposes  of  comparison.  This 
is  the  standard  algorithm  for  determining  the  best  solution  of  an  overdetermined  system  of  linear  equa¬ 
tions  (2).  It  minimizes  the  sum  of  the  squares  of  the  distances  to  the  measured  bearing  lines.  Imple¬ 
menting  this  procedure  does  not  require  that  all  of  the  angle  measurements  be  stored  in  memory  but 
rather  five  sums  are  accumulated  in  memory  as  the  various  fixes  are  taken.  Based  upon  this  informa¬ 
tion  an  estimate  of  emitter  location  is  obtained  by  solving  a  system  of  linear  equations. 

For  the  particular  set  of  conditions  one  is  concerned  with  hers,  the  combined  effect  of  large  angle 
errors  and  short  base  line  causes  the  standard  algorithm  to  sometimes  fail  to  even  approximate  titer 
correct  emitter  position,  whereas  these  new  algorithms  do  work.  The  new  algorithms  approximate  the 
condition  that  minimizes  the  sunt  of  the  squares  of  the  angle  errors.  Although  the  exact  imposition  of 
this  condition  leads  to  numerical  problems  which  are  solvable  only  by  time  consuming  iterative 
processes  t3|,  the  approximation  presented  herein  gives  a  new  algorithm  compilable  to  the  standard 
algorithm  in  that  it  also  requites  the  accumulation  in  memory  of  only  five  sums.  The  estimate  of  the 
emitter  location  is  again  determined  by  a  system  of  linear  equations. 

In  the  course  of  deriving  this  latter  algorithm,  an  intermediate  algorithm,  referred  to  as  the  qua¬ 
dratic  algorithm,  will  be  arrived  at  whereby  emitter  location  is  determined  by  a  simultaneous  system 
consisting  of  a  linear  equation,  a  quadratic  equation,  and  a  sign  condition.  With  further  reckoning,  this 
system  is  reduced  to  a  pair  of  linear  equations  whereby  an  algorithm  referred  to  as  the  asymptotic  algo¬ 
rithm  is  determined  The  last  algorithm  given  is  the  result  of  extending  the  asymptotic  algorithm  so 
that  it  determines  emitter  velocity  as  well  as  position. 

live  general  configuration  being  considered  is  that  of  an  emitter  along  with  a  moving  ESM  plat* 
form  from  which  angle  measurements  on  the  emitter  are  made  at  a  number  of  positions  indexed  by  i 
which  runs  front  1  to  A'.  Notation  used  in  presenting  the  algorithms  referred  to  above  is  introduced  as 
follows  and  is  illustrated  in  Fig.  1: 

4  Perpendicular  distance  from  estimated  emitter  position  to  the  line  of  sight  for  the 
i  th  bearing  measurement 

Position  of  ESM  platform  for  the  i  th  took  at  the  emitter 
(x,y)  Cartesian  coordinates  of  the  estimated  emitter  position 
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Fig.  I  —  Geometric  configuration  of  the  emitter  anti  moving 
platform  illustrating  the  angle  measurement* 


(r.d)  Polar  coordinates  of  the  estimated  emitter  position 

ill ,  Emitter  bearing  as  measured  by  the  i  th  look  from  the  ESM  platform 

fob ,  Angular  difference  between  I  th  measurement  and  calculated  estimate 

of  the  emitter  position 

If  the  distances  involved  are  such  that  the  earth's  curvature  is  significant,  the  above  quantities 
are  defined  in  a  plane  tangent  to  the  earth's  surface  at  a  point  near  the  platform  position.  The  resulting 
emiuer  position  is  then  projected  from  the  plane  onto  the  earth’s  surface. 

THE  DISTANCE  LEAST  SQUARES  ALGORITHM 

In  this  section,  the  algorithm  that  minimizes  the  sum  of  the  squares  of  the  distances  to  the  meas¬ 
ured  bearing  lines  is  derived.  The  distance  to  the  f  th  bearing  line  is  given  by 

4  *“  la,x  +  4y  “  i; I*  0) 

where 

o;  -  $int£  „  (2) 

4** ‘■•eosdi,.  0) 

and 

c(-  p,  sin  ip,  -  qt  cos*,.  (4) 

The  quantity  to  be  minimized  is  expressed  as  follows: 

M  -  £  41.  (5) 

>-! 

Taking  partial  derivatives  of  M  with  respect  to  x  and  y  gives  the  following  system  of  linear  equations  for 
these  variables. 
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B[X  +  Ay  *=  E 
Ax  +  =  -  D 
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(6) 

(7) 


Here,  the  following  notation  has  been  incorporated: 
,v 

A  -  X  a>  b‘- 

/-i 

flt  -  X  a/2. 

/*»  t 

fl»-  £  b,\ 

/-i 

d  «*  -  X 


and 

E-tm. 

(»t 


(8) 

(9) 

(10) 

(U) 

(12) 


(13) 


(14) 


This  derivation  determines  the  distance  least  squares  algorithm  for  calculating  the  emitter  posi¬ 
tion.  It  consists  of  the  following  equations  to  be  evaluated  in  the  order  given:  2,  3,  4,  8,  9, 10, 11, 12, 
13,  and  14. 

THE  QUADRATIC  ALGORITHM 

In  tins  section,  an  algorithm  is  devised  that  for  small  base  line  and  large  angle  errors  gives  consid¬ 
erable  Improvement  in  the  accuracy  of  the  results.  The  algorithm  approximately  minimizes  the  sum  fif 
the  squares  of  the  angle  errors  £>& This  condition  is  expressed  as  follows; 

M-  £<W.  (15) 

*»t 


The  simultaneous  solution  of  Eqs.  6  and  7  is: 
‘  BiE  +  AD 
B&-A* 
and 

-#,0“  AE 

u  «  — . — 

B\  li}-A2 


where 


4$ ,  “  sin'  ’ 


____ _ _ _ _ 

yf(x  -  p4)J”+  (y  -  <t,)} 


(16) 


To  arrive  at  a  practical  algorithm  for  determining  the  emitter  position  associated  with  this  minimi¬ 
zation,  one  must  assume  the  following  approximations:  ihc  USM  platform  positions  are  ail  in  the  vicin¬ 
ity  of  the  origin  and  the  angle  4t>,  can  be  approximated  by 

4*,-—.  (17) 

r 
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With  these  approximations  the  quantity  to  be  minimized  is  given  by 

-£l4T. 


Expressed  in  polar  coordinates  </,  is  given  as 


dt  -  I atr  cos#  +  b,r  sin#  -  c,|.  (19) 

Using  this  expression  in  Eq.  18  gives  the  following: 

N  Ci  Y 

At  »  £  a,  cos  #  +  b,  sin  # - .  (20) 

(-t  r 

To  minimize  At,  partial  derivatives  of  the  above  expression  are  taken  with  respect  to  rand#.  The 
following  expressions  result. 

**■  2  £  a,  cos#  +  bj  sin#  -  —  (21) 

i-l  r  r 

*•  2  V  a,  cos#  +  b,  sin #  **  ~  I-  e(  sin #  +  b,  cos#  j,  (22) 

O#  f«|  1  '  3 

The  partial  derivatives  represented  by  Eqs.  21  and  22  are  set  equal  to  zero.  The  resulting  equations  ere 
multiplied  by  appropriate  factors  and  expressed  in  terms  of  Cartesian  coordinates.  These  steps  result  in 
the  following  system,  a  quadratic  and  a  linear  equation  for  x  and  /. 

.4a3  +  B*y  -  Ay3  +  Ox  +  Ey  »  0  (23) 

U  -  Oy  -  G\  (24) 

&  and  G  are  given  os  follows: 

B-Bi-Bi  (25) 


G  -  £  t,K  (26) 

f-t 

The  coefficients  A,  By,  Bj,  D,  and  £  have  previously  been  defined.  A  variable  u  Is  introduced  with  the 
equation 

Ox  +  By  -  a,  (21) 

which  is  solved  simultaneously  with  Eq.  24  to  give 

<j» 

+  d1 


Eu  -  DG 
E3  D1  ‘ 
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Substituting  Eqs.  28  and  29  into  Eq.  23  gives  the  following  quadratic  equation  for  ic. 


Au2  +  flu+r  =  0.  (30) 

Here,  the  coefficients  are  given  by  the  following  expressions: 

A  -  A  (D2  -  E2)  +  BDE,  (31) 

0  -  4ADEG  +  BG(E 2  -  D2)  +  (E2  +  D2)2.  (32) 

and 

r  «  G2(AE2  -  BDE  -  AD2) .  (33) 


Using  the  quadratic  formula,  u  is  obtained  as 
-n  ±  v'n5  -  4Ar 


a  — 


2A 


(34) 


To  determine  the  ±  sign  in  Eq.  34,  one  must  use  the  following  condition  which  assumes  a  minimum 
value  of  M: 

~r  >  o  (35) 

Carrying  out  the  indicated  partial  differentiations  results  in  the  following  condition: 

fix1  -  4Xxy  -  By2  +  Ex  -  Dy  >  0.  (36) 


Hence  the  ±  sign  in  Eq.  34  is  chosen  to  satisfy  Eq.  36. 

This  derivation  determines  the  quadratic  algorithm  for  calculating  the  emitter  position.  It  consists 
of  the  following  equations  to  bo  evaluated  in  the  order  given:  2,  3,  4,  8,  9,  10,  25,  II,  12,  26,  31,  32, 
33,  34,  28,  29,  and  36. 

THE  ASYMPTOTIC  ALGORITHM 


In  this  section  the  quadratic  algorithm  of  the  previous  section  is  simplified  using  an  asymptotic 
approximation.  Equation  23  is  a  hyperbola  which  may  be  approximated  by  one  of  its  asymptotes.  This 
leads  to  a  new  algorithm,  which,  for  all  practical  purposes,  gives  nearly  identical  results  to  those  of  the 
quadratic  algorithm. 

Determination  of  the  asymptote  is  accomplished  in  a  coordinate  system ,  referred  to  as  the  primed 
system,  for  which  the  axes  have  been  rotated  through  an  angle  w.  In  this  system,  Eqs.  23  and  24,  after 
being  transformed,  retain  the  same  form  with  the  corresponding  variables  and  parameters  denoted  by 
primes.  The  transformed  equations  arc 

A  V3  +  BYy'  -AY2  A-  DY  +  EY  “  0  (37) 

and 

EY  -  DY  “  G‘.  (38) 


The  primed  quantities  are  expressed  in  terms  of  the  unprimed  quantities  as  follows: 

x'  -  x  cos  «  +  y  sin  w  (39) 
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-  x  sin  o)  +  y  cos  w 

(40) 

A  cos2  a)  +  B  sin  <a  cos  w  -  A  sin2  w 

(41) 

-  4,4  sin  w  cos  <u  +  0 (cos2  cu  -  sin2  w) 

(42) 

D  cos  <o  +  E  sin  to 

(43) 

-  D  sin  to  +  E  cos  <o 

(44) 

G 

(45) 

The  primed  coordinate  system  is  chosen  so  that  the  asymptotes  are  parallel  to  the  coordinate  axes.  This 
is  accomplished  by  choosing 

d'-O.  (46) 

Imposing  this  condition  on  Eq.  41  gives 

A  cos2  a  +  B  sin  <d  cos  o>  -  A  sin2  w  -  0  (47) 

Front  this  relation  the  sine  and  cosine  of  the  angle  of  rotation  are  determined  to  be 


COS  <a  «•  -  lot  A, 

(48) 

and 

sin  <u  —  a  $  -  B), 

(49) 

where 

(SO) 

and 

0  **  Jaa'*  +  i)2. 

(SI) 

Substituting  Eq$,  48  and  49  into  Eqs. 

39  through  45  gives  the  following  set  of  equations  representing 

the  coordinate  transformation. 

x'-ab 2Ax  +  ((3  -  B)yl 

(52) 

/-« (-0  “  B)x  +  Uyl 

(S3) 

B’~t 3. 

(54) 

V-al-lAD  +  h'ifi  -  fl)|. 

(55) 

and 

£'-u(-£>(0  -  B)  -  2A£T 

(56) 

Applying  the  condition  of  Eq.  46  to  Eq.  37  gives,  for  the  equation  of  the  hyperbola  in  the  primed  sys¬ 
tem, 

8'xY  +  D'x' +  £'/ ~  0.  (57) 


(SS) 


This  equation  may  be  expressed  as  follows: 
(SV  +  r><<8 y  +  M  -  D'£\ 


6 


NRL  REPORT  8483 


Its  asymptotes  are 

3'x'  +  £'  -  0  (59) 

and 

£y  +  D’-0.  (60) 


The  condition  that  M  be  a  minimum  rather  than  a  maximum  makes  Eq.  36  in  the  unprimed  system 
become  in  the  primed  system 


B'W2  -  y'2)  +  £‘x'  -oy>  0. 


(61) 


This  relation  may  also  be  expressed  as 
12  ,  *•*.  , 

/  + 


E' 


28' \ 


B‘ 

IB' 


jy_ 

28' 


D' 

28' 


>  0. 


(62) 


The  asymptote  represented  by  Eq.  60  satisfies  this  condition  for  values  of  *'  far  from  the  origin  where 
the  asymptote  is  a  good  approximation;  minimum  Af  is  obtained  for  *' outside  the  interval  0). 

In  the  unprimed  system  this  asymptote  is  given  by 

Rx  +  Sy~  T.  (63) 


where 

■  '  (64) 

S  -  2A$.  (65) 

and 

r-  BQ  -  B)  -  2AD,  (66) 


The  estimate  of  the  target  position  is  calculated  by  solving  the  simultaneous  linear  system  consist* 
ing  of  Eqs.  24  and  63.  This  results  in  the  following  expressions  for  the  target  position: 


US  4-  07 
ES  +  DR 


m 


and 


t;r  -  gr 

*~~ES  + DR’ 


(68) 


The  algorithm  resulting  from  this  derivation  consists  of  Eqs.  2,  3.  4,  8*  9, 10,  25,  tl,  12,  26,  51, 
64, 65, 66, 67,  and  68,  evaluated  in  tha*  order. 


The  process  followed  in  applying  this  algorithm  may  be  summarized  as  follows.  For  each  bearing 
measurement  <l>.  it  is  necessary  to  know  the  platform  position  (/?,,<?,).  Five  sums  are  accumulated  tit 
memory:  the  quamtlics  a,  B,  A  £,  and  <#,  which  determine  a  pair  of  linear  equations  to  be  solved  for 
the  estimated  emitter  position. 

SAMPLE  RESULTS  FOR  THE  ASYMPTOTIC  ALGORITHM 

Results  are  shown  in  Table  1  for  the  asymptotic  algorithm  under  a  number  of  assumed  conditions. 
The  true  emitter  range  is  290  distance  units.  The  ESM  platform  moves  normal  to  the  uuc  emitter 
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5 

S 

5 

5 

5 

5 

S 

S 

5 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

SO 

50 

SO 

SO 

so 

so 

so 

50 

SO 

so 

so 

so 

so 

so 

100 

100 

100 


Table  1  —  Sample  Results  for  the  Asymptotic 
Algorithm  (range  =  200  units) 


Base 

units) 

Sigma 

DF 

(deg) 

Sigma 

Range 

(units) 

HUSH 

■ 

■Era 

5 

0.1 

17.6 

.04 

8.81 

10 

0.1 

8.6 

.04 

4.28 

10 

0.2 

17.6 

.09 

8.81 

20 

0.1 

4.2 

.04 

2.12 

20 

0.2 

8.6 

.09 

4.28 

50 

0.1 

1.7 

.05 

0.86 

50 

0.2 

3.4 

.09 

1.72 

50 

O.S 

8.7 

.22 

4.33 

50 

17.9 

.45 

8.93 

5 

i  0.1 

13.0 

.03 

6.50 

10 

0.1 

6.6 

.03 

3.28 

10 

0.2 

13.0 

6.50 

20 

0.1 

3.3 

.03 

1.65 

20 

0.2 

6.6 

3.28 

20 

0.5 

16.3 

.17 

8.13 

SO 

0.1 

1.3 

mmp  ‘HHH 

0.67 

50 

0.2 

2.7 

1.34 

50 

OS 

6.6 

.17 

3-31 

SO 

1.0 

13.1 

.35 

;  6.57 

s 

0.1 

11.1 

.03 

S.54 

10 

o; 

5,5 

.03 

2.77 

10 

0.2 

111 

.05 

S.54 

20 

0.1 

2.8 

.03 

1.39 

20 

0.2 

S.S 

.05 

2.77 

20 

O.S 

13.9 

.13 

6.95 

50 

0.1 

U 

.03 

0.56 

SO 

0.2 

2.2 

,05 

1.12 

SO 

O.S 

5.6 

.13 

2.80 

SO 

1.0 

11.2 

.27 

5.60 

S 

0.1 

6.6 

02 

3.31 

s 

0,2 

13.5 

.03 

6.73 

10 

0.1 

3.3 

.02 

1.65 

to 

0.2 

6.6 

.03 

3  31 

10 

O.S 

17.0 

.08 

8.5 1 

20 

0.1 

16 

.02 

0.82 

20 

0.2 

3.3 

.03 

1.65 

20 

O.S 

8.3 

.OS 

4.16 

20 

1.0 

17.0 

.16 

8.52 

SO 

0.1 

0.7 

.02 

0.33 

SO 

0.2 

1.3 

.03 

0.66 

50 

0.5 

3.3 

.08 

1.66 

SO 

1.0 

6.7 

.16 

3.34 

so 

2.0 

13.6 

.32 

6.81 

s 

0.) 

5.1 

2.56 

s 

0.2 

10.4 

5.21 

10 

0,1 

2.6 

1.28 
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Table  1  —  (Continued) 


Looks 

Base 

(units) 

Sigma 

DF 

(deg) 

Sigma 

Range 

(units) 

100 

10 

0.2 

5.1 

0.02 

2.57 

100 

10 

0.5 

13.2 

0.05 

6.59 

100 

...  20 

0.1 

1.3 

0.01 

0.64 

“100 

20 

0.2 

2.6 

0.02 

1.28 

100 

20 

0.5 

6.4 

0.05 

3.22 

100 

20 

1.0 

13.2 

0.09 

6.60 

100 

50 

0.1 

0.5 

0.01 

0.23 

100 

50 

0.2 

1.0 

0.02 

0.51 

100 

50 

0.5 

2.6 

0.05 

1.29 

100 

50 

1.0 

5.2 

0.09 

2.59 

100 

50 

2.0 

10.5 

0.18 

5.27 

200 

5 

0.1 

3.2 

0.01 

1.62 

200 

5 

0.2 

6.5 

0.02 

3.26 

200 

5 

0.5 

16.8 

0.04 

8.39 

200 

10 

0.1 

1.6 

0.01 

0.81 

200 

10 

0.2 

3.2 

0.02 

1.62 

200 

10 

0.5 

8.2 

0.04 

4.09 

200 

10 

1.0 

16.8 

0.07 

8.39 

200 

20 

0.1 

0.8 

0.01 

0.40 

200 

20 

0.2 

1.6 

0.02 

0.81 

200 

20 

0.5 

4.1 

0.04 

2.03 

200 

20 

1.0 

8.2 

0.07 

4.09 

200 

20 

2.0 

16.8 

0.15 

8.41 

200 

50 

0.1 

0.3 

0.01 

0.15 

200 

50 

0.2 

..  0.6 

0.01 

0.32 

200 

50 

0.5 

1.6 

0.04 

0.81 

200 

50 

1.0 

3.3 

0.07 

1.64 

200 

50 

2.0 

6.6 

0.15 

3.29 

200 

50 

5.0 

17.1 

0.37 

8.S5 

500 

5 

0.1 

2.1 

0.01 

1.05 

500 

5 

0.2 

4.2 

0.01 

2.10 

500 

5 

0.5 

10.5 

0.02 

5.26 

500 

10 

0.1 

1.1 

0.01 

0.53 

500 

10 

0.2 

2.1 

0.01 

1.05 

500 

10 

0.5 

5.3 

0.02 

2.63 

500 

10 

1.0 

10.5 

0.04 

5.26 

500 

20 

0.1 

0,5 

0.00 

0.26 

500 

20 

0,2 

1.1 

0.01 

0.53 

500 

20 

0.5 

2.6 

0.02 

1.31 

500 

20 

1.0 

5.3 

0,04 

2.63 

500 

20 

2.0 

10.6 

0.09 

5.28 

500 

50 

0.1 

0.2 

0.01 

0.10 

500 

50 

0.2 

0.4 

0.01 

0.21 

500 

50 

0.5 

1.1 

0.02 

0.54 

500 

50 

1.0 

2.1 

0.04 

1.06 

500 

50 

1  Id 

0.09 

2.12 

500 

50 

5.0 

Mi 

0.22 

5.36 
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direction  along  base  lines  of  5,  10,  20,  and  50  units.  As  the  base  line  is  traversed  the  number  of  equally 
spaced  angular  measurements  is  5,  10,  20,  50,  100,  200,  and  500.  A  Gaussian  random  error  is  intro¬ 
duced  into  the  bearing  i /»,-  for  each  look.  The  standard  deviations  for  these  errors  "Sigma  DP  are  0.1°, 
0.2®,  0.5°,  1°,  2°,  and  5°.  For  each  combination  of  the  above  parameters,  100  Monte  Carlo  runs  were 
made.  Standard  deviations  (over  the  100  runs)  were  calculated  for  the  estimated  range  and  bearing. 
The  average  percent  error  was  calculated  from  the  ratio  of  range  standard  deviation  to  true  range. 
Table  1  shows  results  of  combinations  for  which  the  average  percent  error  was  less  than  10%. 

Consider  the  ellipse  with  semimajor  axis  given  by  the  range  standard  deviation  and  semiminor 
axis  given  by  the  product  of  the  range  and  the  tangent  of  the  bearing  standard  deviation.  The  probabil¬ 
ity  of  the  true  emitter  position  falling  within  this  elipse  is  0.393.  This  probability  becomes  one-half 
when  the  linear  dimensions  of  the  elipse  are  multiplied  by  the  factor  1.177.  For  a  circular  error  proba¬ 
bility  (CEP)  of  one-half,  the  radius  is  obtained  by  multiplying  the  range  standard  deviation  by  the  fac¬ 
tor  0.675. 

A  paper  by  Wegner  [4]  reports  on  the  location  accuracies  that  are  possible  using  optimum  estima¬ 
tion  procedures.  The  errors  indicated  in  Table  1  are  in  reasonable  agreement  with  errors  determined 
under  like  conditions  from  Wegner’s  paper,  suggesting  that  the  asymptotic  algorithm  is  nearly  optimal. 

AN  ALGORITHM  FOR  A  MOVING  TARGET 

An  algorithm  that  obtains  emitter  velocity  as  well  as  position  is  stated  here.  The  algorithm  will 
not  yield  an  answet  when  measurements  arc  made  from  a  platform  in  straight  line  motion.  The  plat¬ 
form  must  make  a  nonlinear  maneuver.  For  this  algorithm  in  addition  to  obtaining  at  each  look  the 
angular  fix  4>i  and  the  platform  position  (ft,?,),  the  time  it  is  also  recorded  and  utilized.  Equations  69 


through  74  represent  the  calculations  for  a  single  look. 

Ai  **  —  sin  i/t  i  cos  i  (69) 

J9„~sinj0,  (70) 

!i2t  *  COS2  »/<;  (71) 

D,  -  p,  sin  <l> ,  cos  ,  -  q,  cq$j  $ ,  (72) 

E,  m  pi  $inJ  >{i ,  -  q,  sin  ,  cos  (73) 

Gj m  (p,  sin  i -  <f,  cos  ,)J  (74) 


Implementation  of  the  algorithm  requires  the  accumulation  in  memory  of  the  19  sums  given  by  Eqs.  75 


through  90, 

£cos»<l»/  (75) 

cl 

Q  “  £$int5,  (76) 

r-l 

4  ~  £  A,  (77) 

i*\ 

A’-J't.A,  (78) 

4,/«  £  (79) 

M 

^*“£4.  *-1.2  (80) 


10 


NRL  REPORT  8483 


BL-it,Bkll  k~  1,2 

<-i 

B?-  I  /,’*«.  k-  1,2 

N 

D-  I  A 

/-I 

£>'-  ff,  A 

/-i 

/>"-  £  tin, 

i- 1 
N 

E-lE, 

/-i 

E'-tt.E, 

/-i 

£  */£, 

/•  l 

c-  Xft 

-f.. 

/*  i 

When  a  position  and  velocity  arc  desired,  the  next 
91  through  98. 

/<  -  fit  <?  ~  /!/> 

R'«  fl(  G- 

£{'<?-  ^"P 
S  —  y#C?  “  Sj  P 
S'-  4'Q  ~  BiP 

Sitm  aHq  „  £^/  p 

r-  EQ  +  DP 
T'~  E'Q  +  D’P 


(81) 

(82) 

(83) 

(84) 

(85) 

(86) 

(87) 

(88) 

(89) 

(90) 

of  formulas  to  be  implemented  is  given  by  Eqs. 

(91) 

(92) 

(93) 

(94) 

(95) 

(96) 

(97) 

(98) 


The  initial  position  (x.y)  and  the  velocity  (v,w)  arc  determined  by  solving  the  following  linear  system. 


Ex  —  Dy  -V  E*v  ~  iPw  «•  G  (99) 

-  D'y  +  £'"v  -  Z)"w  «  G‘  (100) 

Hx  4 •  Sy  +  Rfv  ±  S'w  -  f  (101) 

-f  +  P,V  +  •“  7y  (10?) 


SAMPLE  RESULTS  TOR  THE  MOVING  TARGET  ALGORITHM 

Listed  in  Tabic  2  is  u  number  of  sample  calculations  using  the  moving  turget  algorithm.  The 
emitter  range,  lengths  and  orientation  of  the  base  iincs,  errors  in  the  angular  measurements,  and  the 
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Table  II  —  Sample  Results  for  the  Moving  Emitter 
Algorithm  (range  =  200  units) 


Looks 

Base 

(Units) 

Sigma 

DF 

(deg) 

Sigma 

Range 

(units) 

Sigma 

Bear. 

(deg) 

%  Error 
Range 

%  Error 
Total 
Displ. 

%  Error 
Radial 
Displ. 

%  Error 
Angular 
Displ. 

20 

50 

.10 

4.0 

.12 

2.01 

10.28 

18.24 

2.38 

50 

SO 

.10 

2.9 

.08 

1.45 

7.10 

12.69 

1.61 

50 

50 

.20 

5.8 

.16 

2.91 

14.12 

25.45 

3.23 

100 

20 

.10 

6.5 

.09 

3.27 

10.12 

20  32 

1.25 

100 

50 

.10 

1.8 

.05 

0.92 

4.45 

7.97 

1.12 

100 

50 

.20 

3.7 

.11 

1.83 

8.76 

15.95 

2.25 

200 

20 

.10 

5.1 

.06 

2.57 

7.92 

15.88 

0.83 

200 

50 

.10 

1.4 

.04 

0.7? 

3.43 

6.23 

0.75 

200 

50 

.20 

2.9 

.07 

1.45 

6.78 

12.48 

1.50 

500 

10 

.10 

9.6 

.07 

4.78 

9.19 

19.32 

1.00 

500 

20 

.10 

3.0 

.04 

1.50 

4.91 

9.39 

0.50 

500 

20 

.20 

6.0 

.08 

3.00 

9.77 

IS. 76 

1.00 

500 

50 

.10 

0.8 

.02 

0.42 

2.08 

3.72 

0.47 

500 

SO 

.20 

1.7 

.05 

0.84 

4.16 

7.44 

0.92 

5’2't 

L.  5?  . 

.50 

4.2 

.12 

2.11 

10.37 

18.59 

2.32 

number  of  runs  arc  the  same  as  those  given  in  the  section  ealicci  “Samplo  Results  for  the  Asymptotic 
'Hgot.ihm"  Since  a  nonlinear  maneuver  is  necessary,  the  platform  is  assumed  in  this  case  to  traverse 
the  base  line  first  in  one  utrectior.  and  ,hen  in  the  other,  In  the  same  length  of  time  that  it  takes  to 
make  th1,  maneuver,  the  emit'cr  moving  with  a  uniform  velocity  undergoes  a  displacement  of  50  dis¬ 
tance  units.  This  motion  makes  a  45°  ang'e  with  the  initial  emitter  position  vector.  Results  are  present¬ 
ed  in  Table  2  lot  those  cases  where  the  percent  error  in  the  radial  displacement  is  Jess  than  30%.  In  ad¬ 
dition  to  the  quantities  described  in  the  section  called "Sample  Results  for  the  Asymptotic  Algorithm." 
Table  l  gives  On  percent  errors  for  the  total,  radial,  and  angular  emitter  displacement  as  determined  by 
the  algorithm. 

SUMMARY 

The  problem  of  determining  the  position  and,  in  some  cases  the  velocity,  of  an  nattier  from  a 
moving  platform  equipped  with  pissivc  direction  finding  equipment  has  been  addressed.  From  a  series 
of  known  platform  positions  spanned  by  a  relatively  small  brso  line,  a  number  of  bearing  measurements 
arc  made.  When  the  angle  errors  are  large,  locating  an  emitter  Hy  geometric  triangulation  or  the  dis¬ 
tance  least  squares  method  may  not  gi*  4  acceptable  accuracy,  if  so,  the  results  may  be  brought  up  to 
an  acceptable  level  of  accuracy  by  using  the  asymptotic  algorithm  derived  in  this  rcfvrt  For  a  system 
with  good  direction  fbvdfng  accuracy  and  a  apabhity  for  recording  the  time  of  a  meat  tfemeni,  an  algo¬ 
rithm  has  been  presented  for  determining  both  the  position  and  the  velocity  of  the  e  miter.  This  algo¬ 
rithm  requires  a  nonlinear  mancu--  ?r  by  the  platform.  The  methods  given  here  nave  the  virtue  of 
'•squiring  only  a  modest  amount  of  memory  and  computational  logic,  Typical  results  have  been 
presented  in  tables  from  which  it  is  possible  to  judyc  the  merits  of  the  method  for  a  particular  applica¬ 
tion. 
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